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We study the effect of parametric excitations on systems of confined ultracold Bose atoms in 
periodically modulated optical lattices. In the regime where Mott insulating and Superfluid domains 
coexist, we show that the dependence of the energy absorbed by the system on the frequency of 
the modulation serves as an experimental probe for the existence and properties of Mott insulating- 
Superfluid domains. 
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The tunability of ultracold atom gases in optical lat- 
tices makes them ideal candidates for the investigation of 
quantum many-particle phenomena in different regimes. 
As a paradigmatic example, within the Bosc-Hubbard 
(BH) model, bosonic atoms in optical lattices have been 
predicted to undergo a quantum phase transition from a 
compressible, phase-coherent superfluid (SF) phase to an 
incompressible incoherent Mott insulator (MI) induced 
by increasing the ratio of the on-site repulsion U to the 
hopping kinetic energy J [lj. For an infinite uniform 
system on a lattice the phase diagram in the (J/U,/j,/U) 
plane is composed of lobe-like MI regions with integer site 
occupation in an otherwise SF phase p], Q ■ By control- 
ling the ratio J/U via the intensity of the optical lattice, 
the MI-SF transition has been experimentally observed 
(HQ together with an energy gap of the expected magni- 
tude in the MI phase. Experiments are always performed 
in finite inhomogeneous systems with a fixed total num- 
ber of atoms subject to an external confinement which 
has been numerically established to induce the coexis- 
tence of MI and SF domains HE 00. 

The formation of MI shells of different site occupation 
have been recently detected jj| but still, an experimen- 
tal probe of the co-existence of SF and MI domains and 
of their properties is presently missing. It has been pro- 
posed jlfj that the presence of MI-SF domains modifies 
the low-energy excitation spectrum of the system, finally 
affecting its thermodynamic properties. 

In this paper we analyze the energy absorption due to 
a periodic time modulation of the lattice intensity ex- 
ploiting the parametric resonances of the modes of the 
system. We find this response to be a probe of the MI- 
SF domains, particularly at frequencies that correspond 
to energies lower than the MI gap. For these frequencies 
the energy absorption is only due to excitations of the 
SF regions, typically ring-like domains, and is strongly 
peaked at frequencies related to the transverse size of 
the SF rings. An experimental observation of these peaks 
may then provide evidence for the co-existence of MI-SF 
domains and give an estimate for their width. 



A parametric excitation is a process in which the pop- 
ulation of certain modes is exponentially enhanced as a 
result of a periodic temporal modulation of a parameter 
in the Hamiltonian. Parametric excitations have been 
studied theoretically in harmonically trapped BEC 01 
and for condensates in optical lattices in different spa- 
tial dimensions |l2j | . Very recently a periodic variation of 
the optical lattice potential depth was used to probe the 
response of a ID condensate across the MI-SF quantum 
phase transition 0, Here, the energy absorption 

was measured versus the optical lattice modulation fre- 
quency f2, showing a broad spectrum in the SF phase and 
peaked features in the gapped MI one. Theoretical works 
that followed analyzed the parametric instabilities of Bo- 
goliubov modes |l5l Il6j. In parallel, very recent works 
considered the linear response energy absorption of ID 
systems within the bosonization picture |l7| or with nu- 
merical diagonalization of realistic systems [lil flfll l2fl | . 
In the presence of confinement, the low-energy form of 
the dynamic structure factor has been argued to yield 
signatures of the presence of SF-MI domains [20| . 

The BH Hamiltonian describing a system of bosons in 
a periodic optical lattice is 



H = -J J2{ b l b j+ h - c 
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with b i the annihilation operator of a boson at site i. We 
assume nearest- neighbour hopping of strength Jo, and 
an on-site repulsion Uq. The values of Jo and Uo are 
given by J = 4/ 'y/Z(V / 'E R )^exp[-4y/W/E^]E R and 
U = ir 2 a /3l(V /E R ) 3 / 4 E R [2l[, where V a and I are the 
optical lattice depth and spacing, E R is the recoil energy 
and ao the scattering length. The chemical potential fi 
fixes the particle number and Vi is the potential bias at 
site i due to the external confinement. 

We consider a two dimensional (2D) atom cloud sub- 
jected to an optical lattice. In the 2D plane the confining 
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potential Vi grows from the center of the cloud towards 
its boundaries. For Jq/Uq larger than the critical value 
inducing the MI-SF transition, the entire 2D system is 
a SF while below this critical ratio concentric alternated 
2D ring- like MI and SF domains are formed Q. 

In a compressible SF ring the average occupation per 
site n is in between the integer occupations of the neigh- 
boring MI domains. Assuming the length of the ring is 
much larger than its width L, we can view it as a SF 
stripe closed with periodic boundaries and consider mo- 
mentum eigenstates along the stripe. For weak depiction, 
in momentum space the BH Hamiltonian |JTJ in the SF 
stripe reduces, at mean field |13|j to 



where 7max = (C — B)e k nUo/2E k . To lowest order in A 
we may replace E k by its time independent part E k . 

When the arguments of the exponents get large, the 
fast oscillations effectively kill the time-dependence of a k . 
In the regime fit ^ 1 this is almost always the case, ex- 
cept on the slowly oscillating modes with | fl — 2E k 1 1 <C 1 
which fulfill the single second order equation 



i fl-2E, 







with solution a k = rj k e luJ + + r] k e ZUJ - , where 
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Here ij k are set by the initial conditions a k (t = 0) = 
c k (t = 0) and d k (t = 0) = 7 max c_ k (t = 0) and result in 

coherent superpositions of c k (t = 0) and c_ k (t = 0). In a 
Within the tight-binding approximation narrQW window of energies \fl/2 - £g| < | 7max | the fre- 

quencies w± gain an imaginary part and the occupation 
of the ?7 k modes grows exponentially with time (reso- 
nant modes), corresponding to the resonant excitation 
of pairs of Bogoliubov modes with opposite momenta, 
whose total energy is approximately equal to the exter- 
nal modulation frequency. Outside this window, the time 
evolution of the occupation is periodic. The occupation 
of the Bogoliubov modes, (c k (t)c k (t)) , is then 



where the ground state energy was subtracted and a k 

(valid for the deep lattices implied in the MI domain 
formation) we have, at small k, e k = Jo^ 2 |k| 2 . 
The Hamiltonian J2J) is diagonalizcd via the Bogoliubov 



-k — "'k'-'k "T" u k a -w VV1 
\ [( e k + nU o) l E i + !] and results in 
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k c k c k 
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where c k creates a Bogoliubov quasi-particle of 2D quasi- 

1 /2 

momentum k and energy E k = [(e k ) 2 + 2n[/oe k ] 

To the system described by (0 we introduce a periodic 
modulation of the lattice depth V of the form V (t) = 
V (1 + A sin fit), with A <C 1. Thus, J and U become 
time dependent, and to lowest order in .4 we have e k (t) = 
e k (l + B sin fit) and U(t) = U Q (l + C sin fit), where B 
and C are calculated from the relations below Eq. QJ. 
The resulting Hamiltonian is 
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sinh 2 ( 7 t) [2n k (T) + l]+n k (T) 



sin 2 ( 7 t) [2n k (T) + 1] + n k (T) (9) 



H 



E 

k^0 



-k c -k 



v k {t) 



where (t) 



= K 

+ nUn) H 



h.c, 



(4) 



+ AEk(t), with AE k (t) = 
CnUo] sin fit, and 



B)sin fit 



(5) 



The Hamiltonian Q is no longer diagonal, and a time- 
dependent term creating pairs of excitations with oppo- 
site momenta is introduced, responsible for the paramet- 
ric resonance. Following the method of Tozzo et al., we 
introduce the operators a k = c k exp [i J dt'E k (t')] and 
the Heisenberg equations on 10} yield (we set Ti = 1) 



in the regimes 4 7 2 lax > (fl - 2££) 2 and 4 7 ^ ax < (fl - 
2El) 2 respectively, where 7 = [| 7 max - ~ ^Kf\] ^ 
and n k (T) = [e E k/ k B T _ j s the Bose distribution. 
This result is correct as long as the occupation of the 
excited modes (the condensate depletion) is small such 
that (0) and Q hold. The exponential growth of the 
resonant-modes occupation is compatible with the result 
of Tozzo et al. 0- 

In order to calculate the energy absorbed by the sys- 
tem due to the time dependent modulation of the lattice 
depth, we consider a modulation that is turned on at 
t = and switched off at time t = r, bringing back the 
Hamiltonian of the system to The energy absorbed 
by the system during the time < t < t is 



£(T) = E£k[(4(^ k (r)} 



k#0 



(c k (0)c k (0)> 



(10) 



n 



a k = 7max "-k 
-k = 7max a k f 



i(nt-2 J E k dt) _ g i(nt+2 J E k dt)~ 
i(Ot-2 J E k dt) _ e ~i(nt+2 J E k dt)~ ^ 



where we perform a thermal average with the unper- 
turbed Hamiltonian H . The absorption is then obtained 
by inserting @ into (|10|) , with the summation over k lim- 
ited to the range of energies where oscillations are slow. 
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The absorption is strongly affected by the density of 
states (DOS) of the system, which is crucially dependent 
on the shell structure. Indeed, if the width of the ring L 
is much smaller than its length L xl transverse quantiza- 
tion of the 2D momenta yields a spectrum composed of 
separate ID branches, labelled by the integers m. The 



Jh 



branch has a dispersion 



E k.; 



k 2 + (^) 



1/2 



(11) 



where k is the ID wavenumber along the ring and we con- 
sidered the low-energy phononic dispersion with sound 

1 /2 

velocity c = (2fiUoJ l 2 ) . The DOS is character- 
ized by Van-Hove singularities near k — of the form 
u(E) ~E/ [E 2 - El m ] 1/2 where E g , n = E° k=0 m . In ad- 
dition, in the presence of the SF rings, ID surface modes 
exist at the interfaces between domains, with a linear 
dispersion at low momenta flpj | . Their DOS is however 
smooth and shows no singularity, which makes their de- 
tection difficult with the technique proposed here. 

We now examine the dependence of the energy absorp- 
tion llOfl on the modulation frequency fl. The typical 
parameters used in recent experiments ^| are t ~ 30 
ms and f2 < 8 kHz such that fir S> 1 and the occupa- 
tion of the modes is determined by Eq. JJJ). At a given 
fi, modes with energies in a "resonant window" of width 
27 max around D./2 are exponentially populated in time. 
Thus, by varying Q and considering the energy absorbed 
by the system we can get information about its spec- 
trum. For example, in the presence of SF rings, chang- 
ing f2 in the vicinity of 2E g 1 the resonant window moves 
across the van Hove singularity and the absorption shows 
a peak. The absorption is then suppressed as O is further 
increased past the singularity. A similar enhancement is 
predicted every time the modulation frequency meets the 
condition 2(£ - hwJ) < fi < 2{E gm + | 7max |) for 
some m. Fig. [Jshows the energy absorbed by the system 
versus f2 for a modulation of amplitude A ~ 0.001 in the 
lattice intensity V . For a typical Jo « Uq/45 w 0.4 KHz 
(characteristic to the n = 2 MI lobe) this corresponds 
to the amplitudes B = 0.01 in the hopping strength and 
C ~ 0.001 in the on-site repulsion. 

The plot in Fig. ^ is the main result of this paper; 
a measurement of the energy absorption as function of 
the modulation frequency supplies a direct evidence for 
the existence of narrow SF shells, as can be seen from 
the difference between the dashed and the solid lines. 
We point out that the enhanced absorption shown in the 
latter is a result of the divergence in the DOS together 
with the fact that the exponent 7 max does not vanish at 
k = 0. In the low-energy regime 2nUo m we get 

7max ~ (C - B)E Q km /A which is finite for k = if m + 0. 

The energy difference between the transverse branches 
resulting from one SF ring, AE = £ 9 , m -£ g , m -i ~ irc/L 
is inversely proportional to its width L. Therefore the 




n/u 



FIG. 1: The total energy absorption after a modulation of the 
lattice intensity V with A = 0.001 during 10 msec, plotted 
as function of the modulation frequency Q/Uo for the mixed 
MI-SF regime at Uo/Jo = 45 (solid line) for the inner ring, 
of average density n = 1.5 particles per site, and for the 2D 
SF phase at Uo/Jo = 20 (dashed line) and average density of 
n = 1.5. We considered £=257 nm and L = 101. 



spacing between adjacent absorption peaks in Fig. ^is a 
signature of the ring-like SF domains, and their typical 
width may be estimated from it, given Jo, Uq and /. 

This result may be compared with the expected energy 
absorption in the absence of SF rings. In case Jo/Uo is 
very small, the width of the SF rings gets smaller than 
one lattice spacing, and the SF part of the system ef- 
fectively disappears. Then, the only modes that can be 
excited are above the MI gap and no signature appears 
in low energy absorption. In the opposite regime, when 
Jo/Uo is larger than the critical ratio for the MI-SF tran- 
sition, the entire 2D system is SF with a smooth DOS. 
The dashed line in Fig. ^ shows the energy absorption 
as function of the modulation frequency in this regime. 

For a 3D cloud of atoms the arguments above hold 
in essentially the same form. Again, the presence of 
SF shells with thickness L yield several 2D Bogoliubov 
branches with interband separation proportional to 1/L. 
In this case, the 2D nature of the modes does not lead to 
van Hove singularities but rather to a staircase profile in 
the DOS. As a result, the absorption vs. f2 will not show 
peaks, as in Fig. ^ but steps, from the existence and 
spacing of which one can infer the presence of domains 
and their size. Again, in contrast, a pure 3D supcrfluid 
phase would show a smooth absorption spectrum. 

We now analyze the consistency of our results and their 
applicability to experiments. We choose the parameter 
values to be in a regime close to that described in 
The recoil energy for Sodium atoms in an optical lattice 
of wavelength 21 = 514 nm is Er = 2nh x 32kHz. For 
the generation of two MI rings we need Vb — 15 E R , 
yielding Uq ^ 0.4^ and J ~ 0.01 E R , The single band 
BH approximation holds, since \Uon{n— 1) -C hv, where 
Tiv = yjAE R V is the excitation energy to the first excited 
Bloch band. A typical ring of width L ~ 10 1 and n = 1.5 
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results in a transverse gap of E 1 ~ 0.1 Uq, allowing to 
probe about 10 transverse branches below the MI gap 
Uq. The visibility of the absorption peaks is strongly 
related to the duration of the parametric modulation. 
An estimate for the minimal duration required for good 
resolution of the peaks is t ~ 5^7majc(-^?,i) ~ 10 msec, 
feasible in current experimental setups. 

The enhancement in the population of modes at zero 
temperature is triggered by the presence of the zero point 
fluctuations, that supply a "quantum seed". At finite 
temperature these arc augmented by thermal fluctua- 
tions. A second important fact enabling this amplifi- 
cation is the bosonic character of the atoms. If we con- 
sider the "fermionic version" of our hamiltonian (J2J, i.e. 
the BCS mean-field one, and proceed with the periodic 
modulation, no exponential growth of the occupation is 
expected due to the fermionic nature of the excitations. 

The sharp peaks in Fig. ^ resu lt from the confinement 
of the superfluid to narrow and long strips and from the 
long-living excitation modes that characterize superflu- 
ids. To separate the contribution of these two sources, 
we consider the effect of viscosity on the parametric res- 
onances in a narrow fluid strip. We concentrate on the 
parametric excitation of sound waves when the sound 
velocity is modulated through a periodic variation of the 
effective mass M(t) = Mo(l + B sin fit) and the interac- 
tion constant U(t) = Uq(1 + C sin fit). The equation of 
motion for the density fluctuation is 



h + C,k{t)f>k + u1{t)pk = 



(12) 



where u 2 {t) = c s (t)k, c s (t) cx yJU{t)/M(t), ((t) = 
M(t)/M(t) + T]k 2 , and rj is the fluid viscosity. The 
transformation pk(t) = yk(f) exp [— J* £k(t')/2dt'] maps 
Eq. i|12|) onto a parametrically excited oscillator 



~ 2 l + Asinfl(t 



to) Vk = 



(13) 



where A 



[c - (i - n 2 /2u 2 



kfi )B] 2 + ( V k 2 flB /2c4 ) 2 



u k0 — uj k q — yijn, i to lowest order in B and 
C. The resonant y' k s are exponentially enhanced, 
y k (t) ~ expfimaxt] with ^( max ~ \A&k,o- The re- 
sulting dynamics of the sound modes has the form 
p k ~ cxp [(j max - r)k 2 /2)t] exp[-_BsinfK/2]. 

Density fluctuations will grow exponentially in time 
only if the amplitude of the modulation is large enough 
such that 7 mol is real and satisfies j m ax > ilk 2 / 2. Fur- 
thermore, the width of the energy absorption peaks is 
2jmax. Thus, when the amplitude of the modulation is 
large enough for the formation of the peaks, the peaks 
acquire a finite width, which for the m'th peak must 
be larger than rj{irm) 2 / 2L 2 . When the viscosity is large 
enough such that the width of the peaks is larger than 
their spacings, the peaks arc fully smeared. 



In conclusion, we considered an ultracold cloud of Bosc 
atoms subjected to a periodically modulated optical lat- 
tice which induces parametric amplification of its excita- 
tions. The energy absorption as a function of the modu- 
lation frequency may serve as a tool to detect features in 
the DOS of the system, that are induced by a transition 
from a uniform SF phase for weak lattice intensities to 
a mixed MI-SF regime for stronger ones. This enables 
the experimental detection of the coexistence of MI-SF 
domains and the measurement of their typical size. In 
addition, we demonstrated the broadening of absorption 
peaks due to damping in normal fluids which is absent in 
case the domains arc supcrfluids. 
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Note added: During the final stage of preparation of 
this manuscript two experiments appeared testing the 
shell structure with local probes [23|, |2J] . 
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